Excitations are localized and relaxation is hierarchical in glass-forming liquids 
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For several atomistic models of glass formers, at conditions below their glassy dynamics onset 
temperatures, T , we use importance sampling of trajectory space to study the structure, statistics 
and dynamics of excitations responsible for structural relaxation. Excitations are detected in terms 
of persistent particle displacements of length a. At supercooled conditions, for a of the order of 
or smaller than a particle diameter, we find that excitations are associated with correlated particle 
motions that are sparse and localized, occupying a volume with an average radius that is temperature 
independent and no larger than a few particle diameters. We show that the statistics and dynamics 
of these excitations are facilitated and hierarchical. Excitation energy scales grow logarithmically 
with a. Excitations at one point in space facilitate the birth and death of excitations at neighboring 
locations, and space-time excitation structures are microcosms of heterogeneous dynamics at larger 
scales. This nature of dynamics becomes increasingly dominant as temperature T is lowered. We 
show that slowing of dynamics upon decreasing temperature below T is the result of a decreasing 
concentration of excitations and concomitant growing hierarchical length scales, and further that 
the structural relaxation time r follows the parabolic law, log(r/r ) = J 2 (l/T— 1/T ) 2 , for T < T , 
where J, r and T can be predicted quantitatively from dynamics at short time scales. Particle 
motion is facilitated and directional, and we show this becomes more apparent with decreasing T. 
We show that stringlike motion is a natural consequence of facilitated, hierarchical dynamics. 



I. INTRODUCTION 

The behaviors of supercooled glass-forming liquids 
manifest complex correlated particle dynamics. Signa- 
ture behaviors, which appear below an onset crossover 
temperature, T , include super- Arrhcnius growth of re- 
laxation times with lowering temperature, stretched ex- 
ponential time-correlation functions, and transport de- 
coupling. Palmer et al. lj suggest that these behaviors 
follow from hierarchical dynamics in which excitations 
on one scale facilitate dynamics of neighboring excita- 
tions [2J thereby creating excitations on larger scales. 
This idea is encoded by a class of dynamical models, 
so-called kinetically constrained models (KCMs) [3J 0]. 
Results from one such model agree with experimental ob- 



* These authors contributed equally to this work. 

t Corresponding author. E-mail: chandlcr@berkeley.edu 



servations of signature dynamical behaviors [SHI] ■ KCMs 
presume that excitations are localized and free of any 
significant static inter-excitation correlation. All large 
scale effects are due solely to the nature of excitation dy- 
namics, which is facilitated and directional. This paper 
presents our discovery that for several different atomistic 
models, these features are emergent properties of under- 
lying Newtonian dynamics. 



A. Molecular dynamics simulations demonstrate 
facilitated and hierarchical dynamics. 

For each of five atomistic models at several different 
densities and temperatures, we ask the question: How 
does an atom move a distance a between relatively long 
lived neighboring positions? We answer this question 
by augmenting extensive molecular dynamics simulations 
with methods of transition path sampling [8 . In particu- 
lar, molecular reorganization in a supercooled liquid is a 



rare event, and transition path sampling can harvest un- 
biased ensembles of trajectories exhibiting such events. 
By applying this approach, we find several important re- 
sults. 

First, for a smaller than or not much larger than a 
particle diameter, particle displacements that stick to a 
new position for a significant period of time are associ- 
ated with correlated displacements of only a handful of 
neighboring particles. These persistent changes in par- 
ticle position, which serve as indicators of what we call 
"excitations," are closely related to the micro-strings [9] 
found in earlier computer modeling studies of heteroge- 
neous dynamics in glass forming liquids. At supercooled 
conditions, i.e., at temperatures T below the onset tem- 
perature, T , their occurrence is sparse. They arise in 
localized regions of relatively high mobility whose size is 
largely independent of temperature, and whose spatial 
distribution is that of a dilute gas. 

Second, for a given displacement length a, we find that 
the equilibrium concentration of excitations, c a , has a 
Boltzmann temperature dependence 

c a oc exp [- J a (1/T - 1/T )] , T<T , (1) 

where J a grows logarithmically with a, i.e., 

J a -J a ,= 1 J a \u{a/a l ). (2) 

Here, a is the space-filling diameter of a particle, the 
values of 7 and J a are material-dependent, and 7 is of 
order unity. 

Third, we find that observed super- Arrhenius temper- 
ature variation of transport properties follows directly 
from this logarithmic scaling of J a . In particular, as the 
average distance between excitations grows with lower- 
ing T, the energy scale for relaxation also grows with 
lowering T. Specifically, from Eq. pi), we argue that the 
structural relaxation time, r, obeys 



r = T cxp J 2 (1/T - \/T o y 



T<T , 



(3) 



with J = \fiyjdt Ja, where df is the fractal dimension- 
ality of heterogeneous dynamics. The reference time, r OJ 
depends little upon temperature, and it is of the order 
of the time to reorganize a local arrangement of particles 
in the presence of an excitation, which in turn is of the 
order of structural relaxation times of the liquid at tem- 
peratures above T a . We show that Eq. ([3| holds quanti- 
tatively with df/d « 0.9 and 0.8 for physical dimensions 
d = 2 and 3, respectively. 



B. Results demonstrate the perspective of 
kinetically constrained models. 

Equation pi is the temperature dependent form for r 
that successfully collapses disparate transport properties 
of all fragile glass formers [SJ [TU] . The logarithmic scal- 
ing from which it results, Eq. [2j is that of the d — 1 



East model [TT] and higher dimensional generalizations 
0Q2]. Here, by providing a microscopic recipe for de- 
tecting excitations, computing c a and predicting t, we 
validate East-like KCMs as good models for the dynam- 
ics of atomistic glass formers. 

These technical and quantitative advances seem es- 
pecially significant in light of the implied physical pic- 
ture. In particular, we establish herein that the princi- 
pal growing length that governs the slowing of dynamics 
in a structural glass former is the mean distance sepa- 
rating excitations, £ a oc c a ■ Relaxation requires the 
correlated dynamics of neighboring excitations, and the 
greater their separation, the longer it takes to coordi- 
nate their motions. It is a picture that appears signifi- 
cantly different than an often invoked structural perspec- 
tive originating with Adam and Gibbs [THl HI]- That 
alternative imagines an underlying mosaic of ordered do- 
mains, so-called "cooperative rearranging regions." Re- 
laxation follows from reorganizing these regions [I5H17J . 
and dynamics slows because mosaic domains grow. Our 
findings do not preclude a picture based upon growing co- 
operative rearranging regions, but a link between it and 
the localized excitations that we document would need 
to relate £ a to the size of mosaic elements, and whether 
an operational definition of the latter can accomplish this 
task is unclear. 



II. METHODS 

The abundance of results collected in this paper are un- 
precedented by the current standards of this field. These 
results establish the behaviors of many distinct large sys- 
tems, and demonstrate the commonality of these behav- 
iors over a broad range of time scales. Even so, further 
studies could provide additional documentation and per- 
haps refinements of the conclusions we are able to draw 
from the results presented herein. 



A. We study the molecular dynamics of five 
distinct simple liquid mixtures. 

The atomistic models we have simulated are two- 
component mixtures for which particles of type a interact 
with those of type 7 with pair potentials, u a ~(r), where r 
is the separation of the pair. These potentials are shifted 
and truncated Lennard-Jones potentials, 



u al (r) 



U 



(LJ)/ (C) \ « < ,- 



0, 



r > r, 



(c) 

CK7 

(c) 

07 



where w^ LJ ^(r ; cr, e) = 4e[(cr/r) 12 — (er/r) 6 ]. For each, we 
use TV = 10, 000 particles in total, with fN and (1 - /) N 
being the number of A-particles and B-particles, respec- 
tively. The different models are distinguished by different 



choices of length and energy parameters, of mixing frac- 
tion /, and of particle masses wia and tob- For notational 
ease, we use dimensionlcss quantities throughout, where 
Boltzmann's constant is unity, and the unit of time is 
(e/mcr 2 ) -1 ' 2 , where e = e\A, m — mA, and a = ctaa are 
the units of energy, mass, and length, respectively 

The first three models detailed in Table U have been 
developed and studied by others before us. The last two 
are modifications of models examined by Harrowell and 
co-workers [2ll [22]. In particular, we replace the inverse- 
twelve power potential of Refs. [5T| and [22] with Weeks- 
Chandler-Andersen [23 (WCA) repulsive potentials. 

By exploring the dynamics of these models for larger 
system sizes and longer times than have been considered 
before, we find that both variants of the 2D-50:50 sys- 
tem appear to coarsen and freeze. As such, we use those 
systems for illustrative qualitative studies over timescales 
smaller than coarsening times. For quantitative behav- 
ior of reversible transport and related properties, we use 
the other four models, none of which exhibit coarsening 
at the conditions we consider. Our molecular dynamics 
for those systems is carried out with N, V and T fixed, 
using the HOOMD-Blue simulation code [21] on graphics 
processors. 



B. Excitations are manifested by particle 

displacements that persist for significant periods of 

time. 

While glassy systems tend to be locally rigid or 
jammed, some locations are atypically less rigid and pro- 
vide opportunities for structural reorganization. These 
locations coincide with defects or excitations in an other- 
wise jammed material. Their presence can be inferred or 
detected by observing non-trivial particle displacements 
associated with transitions between relatively long-lived 
configurations. Because these displacements are reliable 
recorders of excitations, we often treat them as synony- 
mous. Nevertheless, the two are distinct. Displacements 
refer to dynamics in small parts of trajectories, while 
excitations refer to underlying configurations or micro- 
states. We further comment on this point later in this 
paper. 

Non-trivial displacements - the recorders of excita- 
tions - are different than fleeting vibrations. Inherent 
structures [2FJ US] help distinguish the former from the 
latter. The inherent structure of a configuration of an 
iV-particle system is obtained by steepest descent to the 
nearest minimum in the potential energy landscape |27j . 
The inherent structure evolves as dynamics progresses, 
but most intra-basin vibrations are not apparent in that 
evolution. Excitations are present when the dynamics 
produces significant changes in inherent structure. The 
same effect of removing relatively high-frequency vibra- 
tions can be obtained by coarse graining coordinates over 
a small increment of time. We have generally chosen this 
latter procedure. But we have checked in a few cases 



that our principal results for large enough length scales 
and time scales are unaffected by this choice, provided 
the coarse-graining time is of the order of the typical 
period for atomic length scale intrabasin vibrations and 
this time is of the order of, but generally smaller than, 
the typical instanton time At, defined below. Where rj(£) 
denotes the instantaneous position of the ith particle in 
the system at time t, we use Yi(t) to denote its corre- 
sponding position in the time-coarse-grained structure, 
i.e., 



h(t) = 



61 



dt' r(t + t'), 



(4) 



where St is the coarse-graining time. On the occasions 
where we refer to inherent structure coordinates, we use 
rp (t) to denote the position of particle i in the inherent 
structure at time t. 

Figure [T] renders part of a trajectory from one of the 
models we have considered. It illustrates typical behavior 
of these processed coordinates. While inherent structures 
change discontinuously as the system moves between po- 
tential energy basins, the raw and coarse-grained coor- 
dinates are continuous. The particular trajectory illus- 
trated in Fig. [l|a) shows the appearance of the type of 
motion with which we detect an excitation. It is char- 
acterized by three periods: During the first, the inher- 
ent structure of a tagged particle remains near its ini- 
tial point; following this sojourn, the system undergoes 
a rapid transition that lasts up to a time At, termed 
the "instanton" time, after which the tagged particle is 
found in a distinctly different position for another rela- 
tively quiescent period. 

This behavior suggests using the following functional of 
path as the indicator that at time t particle i is associated 
with an excitation with displacement length a: 

hi(t,t a ;a)= f[ 6(\r i (t + t')-f i (t-t / )\-a) . 

t'=t a /2-At 

(5) 
Here, 0(x) — 1 or for x > or < 0, respectively. The 
products are over time steps of a trajectory that extends 
for a time t a > At. This length of time is at least as long 
as the time for trajectories to commit to one basin or the 
other while traversing a transition state region. In the 
parlance of rare-event dynamics, t a is a plateau or com- 
mitment time [281 129] - a time for which the equilibrium 
average of hi(t,t a ;a) grows linearly with t a . It is typi- 
cally 3 or 4 times the mean instanton time for that dis- 
placement length a, (At) a . As such, the indicator given 
by Eq. [5] discards trajectories that exhibit displacements 
with length scale a if they fail to persist for at least as 
long as the typical instanton time. The instanton time 
is the shortest time separating the initial and final so- 
journs. It is a functional of path, and it varies from one 
trajectory to another. 

Fig. [Tib) illustrates the spatial arrangements of these 
excitations. In particular, at supercooled temperatures, 



TABLE I: Parameters for model systems used in this study." 
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"Dimensionless units employed. See text. 

6 Kob-Andersen model from Ref. 1181 . 

c Wahnstrom model from Ref. |19| . 

d Weeks-Chandler-Andersen repulsive force mixture from Ref. 1201 . 

e WCA potentials with the same particle-size ratios and mixing rules as the inverse power potentials of Ref. |21| 

■^ WCA potentials with the same particle-size ratios and mixing rules as the inverse power potentials of Ref. 1221 



for reasonable choices of a, the pictures give the impres- 
sion that the mean excitation concentration is small. As 
temperature is lowered, typical sizes of excitations are un- 
changed, but their concentration clearly decreases. This 
concentration is recorded by 



N 



— ^2hi(0,t a ;a) 



Vt a 

1 
Vt n 



N 



53e(|r i (t o + <At)„)-f i (0)|-o)),(6) 



where the angle brackets denote equilibrium ensemble 
average, V is the volume and N is the total num- 
ber of particles in the system. This second (approxi- 
mate) equality in Eq. p] where hi(Q, t a ;a) is replaced by 
(|r-t(io + (Ai) a ) — fi(0)| — a), holds provided t a is in- 
deed a plateau time. At conditions where the liquid is 
not supercooled a separation of time scales need not ex- 
ist, and thus there is no plateau behavior. We replace 
/ii-functions with the ^-functions when calculating c a at 
supercooled conditions. 

This quantity c a is the average number of displace- 
ments per unit space-time. We refer to it as the average 
excitation concentration. With transition path sampling, 
without the simplifying approximation of Eq. [6l we col- 
lect satisfactory equilibrium statistics for this average by 
accepting or rejecting trajectories with the weight func- 
tional formed from the indicator functional ft,(i, t a ;a) 
times the equilibrium distribution of trajectories, each 
trajectory having a duration of t a . This is done by first 
equilibrating the system with a trajectory that runs for 
many structural relaxation times. A few short sections 
of this trajectory are then chosen to provide first ex- 
amples of a transition or excitation associated with the 
tagged particle. From these first trajectories, a series of 
shooting and shifting moves are then performed [5J [3D] 
generating an ensemble of thousands of independent ex- 
amples of excitations. The ensemble produced by this 
Monte Carlo walk through trajectory space is the sub- 
space, with proper statistical weight, of those trajecto- 
ries in the equilibrium distribution that exhibit excita- 



tion dynamics [3U]. The commitment time t a is al- 
lowed to vary in transition path sampling [31] so that 
dynamics determines a distribution of these times with- 
out pre-conceived notions of its typical values. In some 
cases, we have used transition path sampling to sample 
excitations in very cold systems out of equilibrium. These 
systems are initialized from configurations of a warmer 
equilibrated system, but with temperatures chosen from 
a Maxwell-Boltzmann distribution at a lower tempera- 
ture. 

Considering the range of possible displacement lengths, 
a should be large enough to record significant displace- 
ments. In addition, values of a should be smaller than 
those for which typical trajectories would likely visit in- 
termediate states for long periods of time. Larger val- 
ues of a would obscure separations in time scales, and 
the equilibrium average (hi) would fail to exhibit linear 
growth with respect to t a [32]. For the systems that we 
consider here, taking a to be 0.2 to 2 particle diameters 
proves to be satisfactory, and for this range of displace- 
ment lengths there is a range of suitable commitment 
times where c a is independent of t a . At supercooled con- 
ditions, the typical mean values and fluctuations of At 
and t a are of the order of 10 2 to 10 3 integration steps, 
and very much smaller than structural relaxation times. 

Transition path sampling [5] |3UJ is a necessary tool in 
this study to extract information at very low temper- 
atures. At moderate supercooling conditions, however, 
where sufficient numbers of excitations are present at any 
one time frame, much of what we have done can be done 
with straightforward molecular dynamics. Indeed, this 
is the procedure we use to evaluate c a from Eq. [6] But 
in doing so, it is necessary to identify the duration of a 
time frame, t a or At. Making that identification with 
transition path sampling is reasonably easy because the 
equilibrium distribution for instanton times can then be 
generated automatically |31j . One can then limit the 
search for excitations with At values not far from most 
probable values. For a given displacement length a, es- 
pecially for a larger than a particle diameter, taking At 
much larger than the most probable value will lead one 
to harvest sequences of several temporally separated ex- 



C\J 



1.5 - 

- 1 ... 



! t 

£ ra 0.5 
f 1 



T 



— raw coordinates 

— coarse-grained 
-- inherent structure 



At 



-4 




T=0.5 



T=1.5 



T=1.7 



T=T =2.1 




>1.0 



0.8 _|_ 
0.6 < 



0.4 § 
< 



0.2 



0.0 



FIG. 1: Excitations can be identified from dynamical processes localized in space and time. Here, identifications are illustrated 
with rendered representative trajectories for the 2D-68:32 system (models in Tabled). The unit of length is a and the unit 
of time is (e/mcr 2 ) 1 ' 2 ; i.e, see text, (a) The upper panel shows the displacement for the unprocessed coordinate of a tagged 
particle at time t together with that displacement in the inherent structure coordinate and in the time coarse-grained structure, 
i.e., |fi(i) — fi(— At/2)|. Here At is the instanton time and t a is the plateau time. The averaging window used is 5t — 0.6. 
The snapshots are drawn to show the space-filling sizes of particles and are colored to indicate the distance of a particle from 
its position at t — —At/2. The arrows indicate the direction of that motion, (b) The lower panels show similarly rendered 
snapshots for the full 10,000-particle system over one instanton time for a tagged particle at the center of the box. Enlargements 
show detailed structures of the excitation dynamics, where arrows indicate the direction of the color-labeled displacements. 
Blue indicates relative immobility while red indicates relative mobility. The depictions are for particle density p — 0.75 at four 
different temperatures (as indicated). The glassy dynamics onset temperature for this system at this density is T « 2.1. 



citations at smaller length scales. With straightforward 
molecular dynamics, one can be assured of not using too 
large a value of At by checking that this time is not be- 
yond plateau value times, as discussed below. 

Many features we highlight in this study of atomic glass 
forming liquids are also features of granular media [33] . 
In the context of granular media, these features have 
been studied with order parameters that focus on dis- 
tinct changes in cages surrounding tagged particles [34] . 
Those order parameters might be useful alternatives to 
those considered in this work. 



III. LOCALIZED EXCITATIONS AND 
HIERARCHICAL DYNAMICS 



With the operational definition of excitation dynamics 
given above, the structure and energetics of this dynam- 
ics can be examined. Thus, in this section we estab- 
lish that the equilibrium statistics of excitation density 
is that of a dilute gas, and that the energetics of ex- 
citations grows logarithmically with the length scale of 
excitation displacement. This logarithmic growth is the 
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FIG. 2: Excitation dynamics exhibits separation of 
time scales with instanton times and plateau times. 

(a) The probability distribution of the instanton time, P a (At), 
for the W system (see Table fil at the lowest equilibrated tem- 
perature for three different displacement lengths, a = a/4, 
a/2 and a. (b) The mean value of the rate function R(t a ) 
demonstrating that a plateau value exists for T smaller than 
the onset temperature, T . The particular system considered 
here is the W model at the density p = 1.296, where T ~ 0.88. 
Lines in Panel b are drawn through data points as guides to 
the eye. 



signature of a particular class of hierarchical dynamics, 
and we show that it predicts the temperature variation 
of structural relaxation times of a supercooled liquid. 



A. Dynamic indicators of excitations are localized 

with spatial and temporal extents that are 

temperature independent. 

The time durations of the dynamical processes that 
produce persistent displacements are characterized by 
the probability distribution of the instanton time, 
P a (At). The behavior of P a (At) for a = a is plotted in 
Fig. Kfa) for temperatures ranging from the onset tem- 
perature to the lowest accessible temperature in our equi- 
librium molecular dynamics simulations. The physical 
meaning of the onset temperature is seen explicitly in 
Fig. [2jb). As T goes below T , molecular motions are 
activated, and the concomitant separation of time scales 
is apparent from plateau values show in that figure. Be- 
yond the plateau region, the function 

R{ta) = p(e(|r i (t a + (At) a )-r i (0)|- a )) 



will decrease with increasing t a . The range of the plateau 
region (i.e., the separation of time scales) increases with 
lowering T. For T close to or above T OJ however, mo- 
tion is not activated and a plateau value does not exist. 
(The determination of onset temperatures is presented 
in the next section). The existence of plateau behavior 
in R(t a ) implies that statistically dominant motions are 
instantonic [28] HI] i like that illustrated in Fig. fTJ 

Transport properties vary by three to four orders of 
magnitude over the range of temperatures considered in 
Fig. Km a), but the distribution of instanton times varies 
little, and its width is far smaller than structural re- 
laxation times at deeply supercooled temperatures (i.e., 
T < T ). In this sense, the extent of excitation dynamics 
in supercooled materials is relatively localized in time. 
Changes in displacement length, shown in Fig.^Ta), pro- 
duce significant changes in P a (At). These changes re- 
late to the hierarchical nature of the dynamics that we 
demonstrate in subsequent sections, where we also report 
on structural relaxation times. 

Along with being temporally localized, excitation dy- 
namics are spatially localized. To show this, we consider 



H(r,t,t';a) = 



(M0,t a ;a)> 



iV 



/ii(0,t a ;a) X)l'*(0-r*(*)l 



ijti 



x<J(fi(i)-fi(t)-r)), (8) 



for t = —t a l1 and t' = t a /2, and with t a in the plateau 
regime. The function pt(r,t,t']a) is the mean displace- 
ment density at r for the time frame t to t' given an 
excitation is at the origin. In the limit of large r, 
j«(r, t, t';a) —> pp^{t — £'), where p is the mean particle 
density and n{t — t') — (|fj(t) — r,(£')|). Deviations from 
this asymptotic limit reflect the degree with which par- 
ticle displacements at r are correlated to the excitation 
dynamics at the origin at time 0. 

This function has some similarity with the ^-functions 
often used to to characterize dynamic heterogeneity (351 — 
157] . Those functions, or susceptibilities, measure mean- 
square fluctuations of densities of particle-position over- 
lap or displacement after a specified time frame. For ex- 
ample, for the special case of taking t = and t' — t = t a , 
the function fj,(r,t,t';a) would be similar to the distinct 
particle contributions to a X4-function at a time t a . But 
the /i-function is different because it is not limited to 
that choice of time variables. It is also different because 
it refers to particle coordinates that are coarse grained 
over time, while X4"f unc ti° ns refer to raw coordinates. 
Further, the /^-function refers to a displacement density 
conditioned on an excitation at the origin, and an exci- 
tation is distinct from a displacement because an excita- 
tion produces a long-lived displacement and not simply 
a fleeting motion. These features that distinguish the pt- 
function from ^-functions are important to the precision 
of our analysis. 
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FIG. 3: Excitations are localized and follow Boltzmann statistics of a dilute gas. (a) Distribution of instanton times for 
excitations with displacement length a — a for the W model at several temperatures below its onset temperature, T ~ 0.88. 
(b) Excess mobility density a distance r from an excitation, relative to its value at r = a for a displacement length a — a. The 
system considered here is the KA model at density p = 1.2 at equilibrium conditions in the supercooled regime, (c) Same as 
(b) but for very cold non-equilibrium conditions, sampled using transition path sampling. The insets of (b) and (c) show the 
corresponding temperature dependence of p(t a ) — (\fi(t a ) — fj(0)|). (d) Temperature dependence of excitation concentration 
c a for the 12 systems considered. The inset shows plots of — lnc a versus 1/T — 1/T for excitations on a length scale a = a. 
The data collapse to the single curve when scaled by the energy parameter J a - Values of T and J a are given in Table [ill 



To use this function as a recorder of excitation size, 
the relevant time frame should surround time and be 
of width t ai so that a pertinent correlation function to 
consider is p{r, — t a /2,t a /2\ a) — pp,(t a ). This function 
oscillates as a function of r manifesting local packing 
of molecules. The same oscillations are found in the 
radial distribution function, (p(r))o , where p(r) is the 
net density in the time coarse-grained structure of parti- 
cles at r, and the average (• • • )o is taken with an ad- 
ditional particle fixed at the origin. Throughout the 
supercooled regime, this radial distribution function is 
essentially temperature-independent, and its correlation 
length (the coarse-graining length over which its oscilla- 
tions disappear) is of the order of one particle diameter. 
Thus, for the purpose of viewing the spatial extent of 
a single excitation, it is useful to divide out the radial 
distribution function and consider 



F(r;a) = 



li(r,-t a /2,t a /2;a) 
(p(r))oM(*o) 



1 



(9) 



The quantity F(r; a) is plotted relative to its value at 
r = a in Fig. p^b). The graphs show that, throughout 



the supercooled regime of the KA model, motions in the 
same time frame correlate over only a small range of dis- 
tances - no more than a few atomic diameters - and this 
correlation range is independent of temperature. This 
observation holds for all equilibrium state points stud- 
ied, Fig. p[b), as well as non-equilibrium state points 
sampled using transition path sampling, Fig. kmc). The 
same holds for all other models we have studied. 

The fact that motions in the same time frame corre- 
late over only a small range of distances implies that the 
underlying configurations, which facilitate motion, must 
also be small with pair correlation lengths no larger than 
the motional correlation lengths seen in Panels (a) and 
(b) of Fig. [3] A large underlying excitation or large corre- 
lation length between such excitations would imply large 
motional correlation lengths, and these are not observed. 
In contrast, as we will see, much larger length scales that 
grow with decreasing temperature are associated with 
structural relaxation and dynamical heterogeneity. In 
this sense, therefore, elementary excitations are spatially 
localized. Growing length scales of structural relaxation 
must arise from correlations between those excitations at 



different time frames. These dynamical correlations are 
hierarchical, as we discuss next. 



B. Excitation densities obey Boltzmann statistics 

with energy scales that grow logarithmically with 

displacement length. 

We have computed the concentration of excitations c a 
from the second equality of Eq. (pf and we find that c a 
obeys the Boltzmann distribution of Eq. (fTl). Figure [3ld) 
illustrates this finding. For the length scale considered 
in that figure, a = a, the plateau time t a is in the range 
20 to 30. The proportionality constant of Eq. (IT]) is of 
the order of (t a a d )~ 1 , changing slightly from one system 
to another. The onset temperature, T , signals the high 
temperature boundary for where c a is characterized by a 
single energy scale J a . We will see that this same tem- 
perature is the high-temperature end to all signatures of 
supercooled glassy dynamics. The energy scale J a de- 
pends upon a in a hierarchical way, Eq. (pi). This finding 
is illustrated in Fig. Ba) . The specific values of J a and 
7 that summarize these results are given in Table [IT) 

We provide two different columns of tabulated data 
for J a because this quantity can be obtained by either of 
two ways - by fitting data to Eq. (n} or by fitting data 
to Eq. pi. It is significant that the two methods give 
similar values. Those obtained in the former way have 
the better statistical certainty, so we use those values for 
predicting relaxation times, which we turn to now. 



C. Relaxation times can be predicted from 
excitation energy scales. 

The logarithmic growth of J a with respect to a has an 
important implication with respect to the temperature 
dependence of transport properties. In particular, and 
in contrast with excitations simply diffusing as a random 
walker, the logarithmic growth implies a hierarchical dy- 
namics as imagined by Palmer et al. [T] in which excita- 
tions on one scale combine to build excitations on larger 
scales. To see how this implies a specific temperature de- 
pendence of transport, consider excitations of displace- 
ment length a for a time slice with thickness of order t a . 
The spatial volume V will be occupied by Vc a jv of these 
excitation displacements, where 1/v is of the order of 
t a a d . The rate at which a given excitation will dissipate, 
1/T a , is the rate at which that excitation can connect 
to neighboring excitations of that same displacement a, 
and the energy to build that connection is its activation 
energy. Therefore, 



l/r a « v exp [- [J la - J a ) (1/T - 1/T )] , T < T 
where 

l a /a=(caa d )-V* 



is the distance to connect a neighboring pair of excita- 
tions. Motions in East-like models are nearly linear [4j[5], 
so that the fractal dimensionality, d{, is expected to be 
close to the physical dimensionality d. It can be less than 
d to the extent that paths connecting excitations are not 
linear. Equations ( 10 ) and dTTj) together with p| yield 



exp 



Jl ( 7 /df) (1/T - 1/T Q ) 2 ,T<T . (12) 



(10) 
(11) 



At a s=s a, this relaxation time is the structural relaxation 
time for the systems we consider. Hence, we find Eq. pL 
with J = J a y/l/dt ■ 

To check this prediction we have determined structural 
times through calculations of the self correlation function 

F s (fe,t) = <exp{ik-[r i (t)-r i (0)]}). 

For the three-dimensional models, we define the struc- 
tural relaxation time, r, according to 1/e = F s (qo,T), 
where go is the wave- vector at which the structure factor 
has its main peak q rj 2w/a. For the two-dimensional 
models, this definition is unsatisfactory because for those 
cases initial cage relaxation makes F s (qo,t) < 1/e long 
before structural relaxation sets in; so for d = 2, we use 
0.1 = F s (q ,T). Figure mb) , shows the excellent agree- 
ment between measured structural relaxation times and 
theoretical prediction for all supercooled systems stud- 
ied. The reference structural relaxation time, r , is de- 
termined by computing a single relaxation time in the 
moderately supercooled regime and using this value to 
align the curves. These values, listed in Table [TTJ are 
approximations to the relaxation times of their respective 
liquids at the onset to supercooled behavior. 

The fractal dimension df can be computed directly |38] 
by clustering neighboring particles that have displaced a 
distance a — a over a time window r. The value of d{ 
is then d{ = dln(n)/dln(r), where n is the number of 
particles contained within a spherical volume of radius 
r positioned at the cluster centroid. We find that this 
procedure gives d{ between 2.5 and 2.7 for physical di- 
mension d = 3, and between 1.8 and 1.9 for d = 2. Any 
values of d{ within those ranges are equally satisfactory 
for fitting the transport data. The values of d{ used in 
Fig. Wb) are df = 2.6 for the three-dimensional systems 
and d{ = 1.8 for the two-dimensional systems, consistent 
with 0.8d and 0.9d, respectively. The one-dimensional 
version of the dynamical scaling we find is that of the 
East model [TT], where d{ = d = I. The remarkable data 
collapse shown in Fig. |4|b) is not the result of simple 
curve fitting. As J a , d{, T and 7 are all determined in- 
dependent of r, the collapse of data is a demonstration 
of the dynamical picture we have derived. 



IV. DYNAMICAL FACILITATION, 

DIRECTIONALITY, AND STRINGLIKE 

MOTION 

The hierarchical nature of structural relaxation that 
we have uncovered in the previous two sections implies 
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FIG. 4: Hierarchical excitation energy scales predict observed structural relaxation times, (a) Demonstration of logarithmic 
growth of J a with respect to a (Eq. [2j) . Fit parameters are given in Table H The inset shows the uncollapsed data. The data 
set 2D-68:32, p = 0.75 is omitted for clarity in the inset, as it scale obscures the other data sets, (b) Structural relaxation 
times predicted from results shown in Fig. [3H and Fig. Eh, which are compared to measured relaxation times. Parameters from 
statistics are listed in Table ITT] Data symbols are the same as those in Fig. [3] 





TABLE II: Parameters for data collapse 


in Figures [3jl andHh. 




Model 


P 


t a 

■I cr 


T 


range 6 


E c 


7 


j d 


To E E 




1.15 


4.0 ±0.3 


0.67 ±0.03 


104 


le-3 


0.55 ±0.05 


4.2 ±0.3 


90 le-2 


KA 


1.2 


5.3±0.2 


0.87 ±0.05 


111 


5e-4 


0.41 ±0.04 


5.7 ±0.3 


120 le-2 




1.25 


7.4 ±0.2 


1.06 ±0.05 


151 


2e-4 


0.39 ±0.03 


7.6 ±0.4 


160 6e-3 




1.3 


8.4 ±0.3 


1.34 ±0.07 


271 


5e-4 


0.39 ±0.03 


9.0 ±0.4 


130 2e-3 




1.2 


4.5 ±0.3 


0.66 ±0.05 


106 


6e-4 


0.59 ±0.1 


4.6 ±0.4 


240 3e-2 


W 


1.25 


6.9 ±0.2 


0.78 ±0.04 


92 


6e-3 


0.42 ±0.03 


6.9 ±0.4 


160 1.2e-2 




1.296 


9.4 ±0.3 


0.88 ± 0.04 


227 


2e-2 


0.36 ± 0.03 


9.3 ±0.3 


180 9.6e-3 


2D-68:32 


0.7 


11.3 ±0.3 


1.36 ±0.05 


51 


3e-3 


0.62 ±0.05 


10.8 ±0.7 


65 6.7e-3 




0.75 


21.8 ±0.3 


2.13 ±0.1 


58 


5e-3 


0.57 ±0.03 


21.8 ±0.5 


60 5e-3 




1.2 


2.1 ±0.1 


0.32 ±0.02 


57 


2e-3 


0.51 ±0.03 


2.2 ±0.2 


60 1.2e-2 


WCA3D 


1.25 


3.5 ±0.2 


0.45 ±0.03 


200 


2e-3 


0.51 ±0.03 


3.6 ±0.2 


135 1.8e-2 




1.296 


5.0 ±0.2 


0.58 ±0.03 


94 


5e-4 


0.47 ±0.03 


5.3 ±0.3 


120 1.8e-2 



"From fits to EqlTlshown in Fig.[3H. 

'The ratio of the maximum to minimum values of c a used in fits 

to Eqll] 
c Error in the fit to Eq [T| given by 1 minus the square of the 

correlation coefficient. 
d From fits to Eqpl shown in Fig. Wk. 
e Error in the fit to Eq [2J given byi minus the square of the 

correlation coefficient. 



that dynamical facilitation [7] plays a central role in 
relaxation phenomena. Dynamical facilitation refers to 
structural rearrangements or excitations allowing for the 
birth and death of excitations nearby in space. The su- 
perposition of such dynamics leads to the formation of 
excitations on longer length and time scales, thus result- 



ing in dynamical heterogeneity [33]. To illustrate these 
features more explicitly, we first consider computer ren- 
dered movies of trajectories, and then turn to quantita- 
tive analysis. 
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A. Facilitation and hierarchical dynamics are 
evident in movies of supercooled liquid dynamics. 

An example of a trajectory is shown in Fig. [5] for 
the 2D-50:50 system. The movie depicts inherent struc- 
tures, with particles colored according to the length of 
their displacement vector over a time window t. Movies 
with time coarse-grained structures are similar, but we 
choose to show those of inherent structures to highlight 
the underlying physics of a most striking feature in these 
movies. These features are the ubiquitous low- frequency 
low-amplitude motions that permeate the system (mani- 
fested by aqua-colored particles). These motions are not 
low frequency harmonic modes, as harmonic motions are 
absent from the inherent structure. Rather, these mo- 
tions are soft, anharmonic motions of the disordered sys- 
tem. On time scales that are long compared to the period 
of these small amplitude soft motions, significant long- 
lived particle displacements occur (yellow, red particles) . 
These motions are associated with what appear to be de- 
fects or excitations in the system. String-like patterns of 
slight mobility (aqua and green particles) surge outward 
from, and retract back towards, the initial excitations [7J. 

Eventually, rare, strong surges cause local regions to ir- 
reversibly deform on larger length scales and time scales. 
These processes, whereby ubiquitous surging excitations 
play a significant role in producing rarer structural re- 
laxations, are an intrinsic feature of East-like models. 
Such surging was noted in Ref. [40] but the connection 
to the East model was overlooked. Activation energies 
grow as displacement length grows, so that frequencies 
of displacements decrease as length scales grow. Further 
larger length scale displacements are built from smaller 
length scale displacements through facilitation. This is 
why in the movie of Fig. [51 it is evident that short time 
local relaxation processes represent a microcosm of re- 
laxation on longer length scales and time scales. As time 
progresses, larger regions of mobility (red, yellow) grow 
outward from the initial relaxed regions, giving rise to 
more collective surges that are characterized by broader 
strings. 



B. Dynamical facilitation, directionality and 

string-like motion are pronounced, and increasingly 

so as temperature is lowered. 

While movies of particle motions are instructive, we 
obtain quantitative measures of hierarchical facilitation 
from the calculation of relaxation time presented in the 
previous section. We can also estimate facilitation vol- 
umes, 



Mt) 



H(r,t a /2,t;a) 
(p(r))o/*(*-W2) 



1 



dr. 



(13) 



The denominator dividing into the /x- function in Eq. ( 13 ) 
is the value of the /x-function in the absence of dynami- 
cal correlations with the initial excitation at the origin. 



We find that the integrand in Eq[l3] looks much like the 
function F(r; a) shown graphically in Fig. p3m and c, but 
with a range that varies with time. Thus, v-p(t) deter- 
mines the volume of space where dynamics at time t is 
correlated to that initial excitation. In the absence of dy- 
namical facilitation, particle displacements are not cor- 
related with the initial excitation at the origin, and v-p (t) 
goes to zero. The qualitative behavior of vp(t) is shown 
in Fig. |6ja) for the W system at several different super- 
cooled temperatures. The other systems we have studied 
behave similarly. 

At low temperatures, the sparsity of excitations makes 
it difficult to obtain accurate statistics for this facilitation 
volume. Irregularities in the variations seen in Fig. K3[b) 
give a sense of the statistical uncertainties in our esti- 
mates for this quantity. The data collected is sufficient to 
reveal the following trends: The facilitation volume ini- 
tially increases as a function of time t, as the initial exci- 
tations facilitate the formation of new excitations nearby. 
vp(t) reaches a maximum value, v™ ax (T), at t ~ r, and 
this maximum value grows with decreasing temperature 
as the distance between the initial excitations increases. 
Over time scales exceeding r, the relaxed regions overlap 
and the system becomes increasingly dynamically uni- 
form with time, causing vp(t) to tend towards zero. Fig- 
ure [61(b) shows that Wp ax (T) consistently grows with de- 
creasing temperature, thus demonstrating the presence 
of dynamical facilitation at all supercooled state points. 

Dynamical facilitation occurs with directionality. This 
is a feature of the d = 1 East model [11] and its d > 2 
counterparts [SJ [12] . For the atomistic models studied 
herein, a quantitative measure of directionality in dy- 
namical facilitation is given by [41] 



N 



u>(r, cos6,t, t ; a) = 

lh J, vr (fti(0,ta;a)^|r«(O-fi(t)l (14) 
(h 1 {0,t a ;a))\ *rr[ 

x <5(cos $i — cos 9) 5{vi (t) — f i (t) — r) \ , 



where 



fj(t,/2)-[fi(t /2)-fi(-t a /2)| 
\fi(t a /2)\ |fi(t„/2)-ri (-t„/2)|' 



(15) 



Displacing particles that lie in the direction of the tagged 
particle's displacement vector have 9i « 0, whereas par- 
ticles that lie in the opposite direction have 64 ~ it. 
A measure of excess oriented motion is therefore |41j 



cv a + ±1 \ u(r,coa6,t,t';a) 
F{r,cos9,t,t ;aj = t-, — —, — — - 1, 



(p(r,9)) of i(t'-t) 



(16) 



where 



N 



(p(r,e)) = {J2d{cos9i-cose)S(f i (t)-f 1 (t)-r) 

Wi / 
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FIG. 5: j jClick to view video clips] ). Time evolution of particle displacements illustrate hierarchical dynamics in the 2D-50:50 
system well below its onset temperature. Specifically, this trajectory runs at a temperature T — 1.1, and the onset temperature 
is T w 2.0. The movie depicts the displacement of each particle in terms of energy-minimized inherent structure coordinates 
\v\ (t) — rj (0)|, with blue indicating overlap with initial position and red indicating a displacement of at least on particle 
diameter. See key. The movie, with frames separated by a reduced time of 25, which is approximately equal to t a for this 
system, spans a time of about 5 r, where the structural relaxation time is r ~ 500£ a • At very early times, non-oscillatory particle 
motions are sparse and spatially decorrelated, indicating the appearance of the initial excitations. As time progresses, regions 
of high mobility (red, yellow) grow outward from the initial excitations until eventually, the regions connect on a timescale 
of about r. Trajectories for the same system at temperatures above the onset temperature do not exhibit these features of 
correlated dynamics. 



12 



2bU 








200 
150 
100 




W, p=1.296, 
a=a, T = 

0.88 0.61 
0.78 0.60 
0.73 0.56 
0.67 0.54 
0.63 0.53 




50 








0.01 


0.1 1 10 100 1000 10000 
t-ta/2 



^ 



E £ 





I 


I ' l 




- 


100 




^Cvf 




- 








,\ 


; 






f^ 




- 


10 




I , I 




; 


I 

( 


) 2 


4 6 




8 



WCA3D, 
p=1.296, 9~n 




J CT (1/T-1/T ) 



4 6 

c (T)-"d,/c (To)- 1/d f 



10 



FIG. 6: Dynamics is facilitated, directional, and string-like, and more apparently so as temperature is lowered, (a) Dynamical 
facilitation volume, vf(£), for a representative W system at several different temperatures, (b) The peak of the facilitation 
volume as a function of temperature for the 12 atomistic models studied, (c) F(r,cos8,t,t';a) for a single temperature with 
parameters a — a, t — t a /2, and t' = 3t a /2 as a function of both the facilitation direction 9 and the distance r from the tagged 
particle at the origin at time t a /2. (d) F(r,cos8,t,t';a) as a function of the temperature with the same parameters as in (c) 
for a single facilitation angle 9 ~ tv. (e) F(r, cos 9, t, t' '; a) for a single temperature with parameters a = a, t = t a /2, and t' — r, 
where r is the structural relaxation time, as a function of both the facilitation direction 9 and the distance r from the tagged 
particle at the origin at time t a /2. (f) F(r,cos9,t,t';a) as a function of temperature with the same parameters as in (e) for a 
single facilitation angle 9 ~ n. (g) Mean string length l s versus the concentration of excitations of displacement length a, c a at 
several different temperatures T for the 10 three-dimensional systems studied. Data symbols are the same as those in Fig. [5] 



Here, as before, (• • • )o denotes the average taken with 
particle 1 exhibiting the excitation at the origin. Fig- 
ure [6tc) shows that on relatively short time scales, exci- 
tations are most likely to occur behind earlier excitations 
a short distance away. Figure Kud) shows that this direc- 
tional preference increases with decreasing temperature. 
Figures W(e) and [6Vf ) show how directional effects dissi- 
pate with increasing time. 

The latter observations are consistent with stringlike 
motion [51 HTM?] . On time scales of the order t a , string- 
like motion takes the form of short "microstrings," whose 
character does not change with the degree of supercool- 
ing [5]. These motions are essentially synonymous with 
elementary excitations that we characterize here. On 



timescales of the order r, microstrings combine to form 
longer strings, with a length scale that grows with the 
distance between the initial excitations. In accordance 
with Refs. [5T], [3] and [17], we construct strings by clus- 
tering particles that fully replace the initial position of 
neighboring particles over a time scale that maximizes 
the mean string length, L s . Strings constructed in this 
way have a minimum size of 2 particles; thus we subtract 
this baseline value from the average when quantifying 
changes in the mean string length, i.e., l s = L s — 2. Fig- 
ure pig) shows that l s grows proportionally as a function 
of the mean distance between initial excitations, c a ■ 
For our comparison, we consider excitations associated 
with particle displacements of length a = a. 
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Earlier ideas about dynamic heterogeneity, specifically 
the concepts of strings EH I43H45] I47H49] and micro- 
strings [5J, seem to be connected with our results. In 
particular, these and other earlier works highlighted the 
presence of correlated stringlikc motions in the heteroge- 
neous dynamics. A small subset of these motions involve 
a few particles that move nearly simultaneously, a pro- 
cess which is called a microstring [9] . Longer strings and 
then clusters of strings are built from these microstrings, 
processes that take place over times that are significantly 
longer than those of a microstring [5] . In the context of 
this paper, microstrings coincide with dynamics of ele- 
mentary excitations, and the building of longer strings 
and clusters of strings coincide with the facilitated hier- 
archical dynamics of excitations. The exponential dis- 
tribution of string lengths reported by earlier studies of 
dynamic heterogeneity in many types of liquids and in 
granular materials [§1 EH EH E3 E3 EH] may now be 
understood as a consequence of the ideal gas statistics 
obeyed by the elementary excitations. Consequently, the 
cooperative length scale of stringlikc motion that grows 
with decreasing temperature coincides with a growing 
distance between excitations. It is a plausible connec- 
tion that is worthy of further study. 



V. DISCUSSION 

This paper shows that excitations are localized, with 
static correlations like those of an ideal gas, and that dy- 
namics in these systems are non-trivial because motions 
occur in a fashion that is facilitated and hierarchical. Per- 
haps most importantly, the number of particles involved 
in a single elementary excitation is independent of tem- 
perature. As a result, dynamics slow and dynamical het- 
erogeneity length scales grow with decreasing tempera- 
ture, not because the fundamental mechanism of particle 
motion becomes increasingly cooperative, but rather be- 
cause the distance between excitations increases. These 
findings would thus seem to conflict with pictures, such 
as the mosaic picture, based on "cooperatively rearrang- 
ing regions" (CRRs) in which the regions themselves are 
considered the fundamental object. Instead, here, we 
view cooperativity as emerging hierarchically over time 
through facilitation, with excitations being the funda- 
mental object. As noted in the Introduction, we do not 
discount the possibility of a yet-to-be discovered link 
that would show how growing CRR size is related to 
the mean distances connecting neighboring excitations, 
4«c a . However, this would presumably require a 
reformulation of the definition of a CRR. The connec- 
tion between string length and mean excitation distance 
demonstrated here might serve as a starting point |50j . 
While we can await that development, at this stage we 
conclude this paper by addressing the most common crit- 
icisms of the picture supported by our results. 

Perhaps the foremost criticism is that no viable proce- 
dure has existed from which localized excitations could 



be seen as emergent properties of atomistic dynam- 
ics [HHHl]. This paper has established such a procedure. 

Another criticism concerns the measured behavior of 
heat capacity. Unlike reversible heat capacity in an equi- 
librium system, this property exhibits significant hystere- 
sis near the glass transition temperature, T g . Heat ca- 
pacity per molecule at temperatures above the range of 
hysteresis is higher than that at temperatures below the 
range of hysteresis by an amount typically larger than 
Boltzmann's constant. In some theories, including Adam 
and Gibbs', this change in heat capacity, AC, is of cen- 
tral importance because it can be related to an imagined 
vanishing of configurational entropy at a finite tempera- 
ture, and this vanishing is supposed to signal loss of er- 
godicity. Inspection of experimental data discredits per- 
ceived correlations between relaxation times and these 
thermodynamic properties (see Fig. 1 of Ref. [10]). Nev- 
ertheless, interest in AC persists, and one may wonder 
whether its typical values are consistent with dynamics 
in structural glass formers pictured in terms of localized 
excitations [5Pj . 

Because localized excitations emerge from molecular 
dynamics, it is unlikely that this picture is inconsistent 
with thermodynamics, at least for the systems studied 
herein. Still, typical values of AC reflect the number of 
translational degrees of freedom removed upon cooling 
below the glass transition, and this number seems phys- 
ically interesting. The removal of translational freedom 
coincides with the removal of fluctuations in the num- 
ber of excitations. Therefore, a simple connection be- 
tween AC and excitation concentration is AC oc c a (T g ). 
This proportional relationship was proposed in Ref. [5], 
where it was shown to be consistent with experimental 
data. The constant of proportionality must grow with 
the number of particles correlated with an excitation of 
unit displacement length. We see from this paper that 
this number is of order 10 or more. Whether further steps 
can be taken to predict precise values of the proportion- 
ality constant remains to be seen. These steps may in- 
volve attempts to understand connections between local 
vibrational modes and excitations. The former are not 
causative of dynamics [52] , but they are correlated to the 
latter |53) , and the former extends over regions of space 
that are large compared to the latter [52] . 

Another criticism is the recent suggestion that facil- 
itation is not the cause of intermittency in dynamical 
activity (births of so-called "avalanches" ) [52 [51] , and 
further that facilitation diminishes as temperature is sig- 
nificantly lowered below (or density is increased above) 
the onset to supercooled behavior [S3]. These seemingly 
contradictory conclusions to our findings are the result of 
differences in the definition of facilitation. In our usage, 
facilitated dynamics is where changes in configurations 
or microstates occur only in the vicinity of excitations. 
Excitations refer to configurations or microstates. Thus, 
in facilitated dynamics, the birth or death of excitations 
occurs only in the presence of neighboring excitations. 
From this definition, it follows that connected lines of 
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excitations permeate trajectory space in systems domi- 
nated by facilitated dynamics [39] . 

Although excitations will connect throughout space- 
time, changes in states will not necessarily connect in 
facilitated dynamics. These changes are the so-called 
"kinks" of kinetically constrained models. Kinks nec- 
essarily record the presence of one or more excitations, 
but excitations can persist for long periods without the 
appearance of kinks. References [32] and [S3] ascribe a 
loss of facilitation to growing spatial and temporal sepa- 
rations between changes in state, but these are growing 
separations between kinks, not excitations. Indeed, the 
behaviors noted in those papers are predicted from the 
facilitated dynamics of the East model and its generaliza- 
tions. Specifically, the motion of particles in the presence 
of excitation lines is intermittent [55) . Clustered bursts 
of activity occur at space-time regions where a particle 
or group of particles intersects excitation lines. Bursts 
end and disconnect from other bursts when this intersec- 
tion ends, after which the particle persists in an inactive 
state for a relatively long period. This behavior is the 
essential element of decoupling phenomena |20j . and it 
becomes more striking as relaxation times grow because 



the lengths of these quiescent periods grow with increased 
supercooling [20, 55J. Whether this understanding can be 
enhanced from analysis presented in Refs. [33] and [S3] 
remains to be seen. 
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